The Single W Production Case 



GlAMPIERO PASSARINO 



Dipartimento di Fisica Teorica, Universita di Torino, Italy 
INFN, Sezione di Torino, Italy 1 



Abstract 

The process e + e~ — > e~V e fif 2 , belonging to the so-called CC20 family, has 
been extensively analyzed in the literature. It is a sensitive probe of anoma- 
lous electromagnetic couplings of the W boson and represents a background to 
searches for new physics beyond the standard model. Moreover, it represents a 
contribution to the e + e~ — * W + W~ total cross section, used to derive a value 
for M w , the W boson mass. The issue of gauge invariance in the CC20 family 
has been solved by the introduction of the Fermion-Loop scheme but several 
subtleties remain, connected with the region of vanishing scattering angle of the 
electron and with the limit of massless final state fermions in a fully extrapo- 
lated setup. A satisfactory solution for computing the total cross section is given 
in the context of the equivalent photon or Wcizsacker- Williams approximation 
which factorizes the flux of quasi-real photons emitted by the electron from the 
interaction rate between the positron and the photon assumed to be real. The 
correct kinematics for the inclusion of initial state QED radiation is established. 
QCD corrections to the process are discussed and numerical results are shown 
and commented. 



e-mail address: giampiero@to.infn.it. 



1 Introduction 



The so-called single W production process is e + e~ — > e~V e W + for small scat- 
tering angles of the outgoing electron and it has been measured at LEP 2 at 
centre-of-mass energies 130 GeV < y/s < 183 GeV using both leptonic and 
hadronic decays of W bosons The signal is, therefore, defined as 

e + e~ — * e~v e ud, (1) 

where stands for a generic up- (down-) quark. In the terminology of it 
is a four-fermion process belonging to the CC20-family. The charge conjugate 
reactions are always understood to be included. The cross section for single W 
production is expected to be small at LEP 2 energies, of the order of 0.5 pb. 
However, this process constitutes a very interesting case both theoretically and 
experimentally. It is a sensitive probe of anomalous electromagnetic couplings 
of the W boson and represents a background to searches for new physics beyond 
the standard model. 




Figure 1 : The CC20 family of diagrams with the explicit component containing 
a fr-channel photon. 

The CC20 process of Fig. [j]is sensitive to the breaking of U(l) gauge invari- 
ance in the collinear limit. For e + e~ — > e~2/ e /i/ 2 , the U(l) gauge invariance 
becomes essential in the region of phase space where the angle between the in- 
coming and outgoing electrons is small, see the work of H and also alternative 
formulations in j|, |) . In this limit the superficial 1/Q 4 divergence of the propa- 
gator structure is reduced to 1/Q 2 by U(l) gauge invariance. In the presence of 
light fermion masses this gives raise to the familiar \sx{m^/ s) large logarithms. 
The correct way of handling CC20 is represented by the so-called Fermion-Loop 
(FL)scheme [[|, the gauge-invariant treatment of the finite-width effects of W 
and Z bosons in LEP2 processes. Briefly, this scheme consists in including 
all fcrmionic one-loop corrections in tree-level amplitudes and re-summing the 
self-energies. However, for practical applications, it has been shown that the 
fixed-width scheme is also satisfactory. Here the cross-section is computed us- 
ing the tree-level amplitude. The massive gauge-boson propagators are given 
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by l/(p 2 + m 2 - iTm). This gives an unphysical width for p 2 > 0, but retains 
U(l) gauge invariance in the CC20 process ||. 

The CC20 process is usually considered in two regimes, |cos#(e~)| > c or 
LACC20 and |cos0(e")| < c or SACC20. Strictly speaking the single W pro- 
duction is defined by those events that satisfy | cos#(e~)| > 0.997 and, therefore 
is a SACC20. 

The LACC20 cross section has been computed by many authors and refer- 
ences can be found in ||. It represents a contribution to the e + e~ — * W + W~ 
total cross section, in turn used to derive a value for M w , the W boson mass. 
This point deserves a comment: by e + e~ — > W + W~ it is meant the ideal cross 
section obtained with the three double-resonant CC03 diagrams and therefore 
the background, FULL - CC03, is evaluated with the help of some MonteCarlo, 
estimating the error on the subtraction by comparing with some other Mon- 
teCarlo. Then M w is derived from a fit to er(CC03) with the help of a third 
calculation. From a theoretical point of view the evaluation of LACC20 is free of 
ambiguity, even in the approximation of massless fermions, as long as a gauge- 
preserving scheme is applied and 9(e~) is not too small. 

For SACC20 instead, one cannot employ the massless approximation any- 
more and this fact makes the calculation unaccessible to most of the Monte- 
Carlos used by the experimental collaborations, with the noticeable exception 
of GRC4F @. EXCALIBUR g is often used in this context with a version 
where a fudge is put so that for one-electron final states on can go down to zero 
scattering angle ^ 

Actually, constructing a CC20 calculation with unconstrained electron scat- 
tering angle is not a problem from the point of view of writing the fully massive 
amplitudes but it is, instead, a question of stability in the numerical integration. 
Moreover the goal of this paper is to show that several subtleties arise in CC20 
for a fully extrapolated setup. 

Single W production and M w measurement are, therefore, complementary. 
Indeed, the phase space requirement | cos#(e~)| > c eliminates events predomi- 
nantly consisting of W pair production since single W production peaks strongly 
at zero scattering angle. 

There is another place where the electron angle cannot be constrained. Ex- 
perimentally events of the type e + e~ — * ud plus a neutrino and an electron, 
possibly in the beam pipe, are not excluded from the hadronic Z lineshape. 
Hadronic events are selected based on final state particle multiplicity in the 
detector, so both genuine high-energy qq events, radiative return events and 
hadronic four-fermion events are selected for the e + e~ — > qq lineshape. This 
gives the total sample. The background is subtracted on MonteCarlo basis, 
using W double-resonant CC03 diagrams, i.e. WW — > all, Z double-resonant 
NC02 diagrams, i.e. ZZ — ► all and hadronic two-photon collisions (e.g. using 
PHOJET g). 

Therefore double-resonant Ws, which are dominant, are treated correctly in 
the experimental procedure and, moreover, single-resonant W's only represent 

1 K. Pittau private communication. 
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a small contribution. The latter could, however, be treated correctly by using 
a CC20 MonteCarlo rather than a CC03 one. In this case the electron is again 
unconstrained and one needs a full angle, massive CC20: the so-called full-CC20, 
or FCC20, regime. 

However, keeping a finite electron mass through the calculation is not enough. 
One of the main results of this paper is to show that there are subtleties in CC20 
also associated with the zero mass limit for the remaining fermions. 

The outline of the paper will be as follows. In Sect. || we introduce the 
general problem of defining the total hadronic cross section at LEP2 energies and 
beyond, and describe how a calculation of FCC20 is to be seen in this context. 
In Sect. H we give a description of FCC20 in the presence of initial state QED 
radiation and show how to implement the correct kinematics for the process. In 
Sect. |] we introduce and discuss the Weizsacker- Williams approximation for a 
small scattering angle of the outgoing electron. The sub-process e + 7 — * V e ud, 
arising in the discussion of the WW-approximation is analyzed in Sects. |^| 
The fully extrapolated setup with massless quarks and QCD corrections are 
presented in Sect. ^. Finally, numerical results and conclusions are shown in 
Sect. 0. 



2 The region of vanishing 9{e ) 

There are at least three applications of CC20 which require an analysis at van- 
ishing scattering angle of the outgoing electron, 9(e~). They are: 

1. the true single W production, i.e. CC20 with |cos#(e~)| > c where, 
usually, c = 0.997, 

2. the evaluation of background for the total e + e~ — > W + W~ cross section, 

3. the evaluation of the inclusive hadronic cross section at LEP2 energies. 

Let us consider in more detail the last application. The FCC20 process is 
not the only background for the total hadronic cross section a(qqX) defined 
as the cross section for qq plus anything. Here, we would like to illustrate the 
general problem, to return in the next section to the study of FCC20. 

The total hadronic cross section that we have defined is an inclusive measure- 
ment of hadron production in e + e _ -annihilation in which production thresholds 
can be seen, e.g. H^-pair production with at least one of the W bosons decaying 
hadronically, or ZZ or other background Jio| . 

Let us repeat what has been done so far in the experimental Collaborations. 
Hadronic events are selected based on final state particle multiplicity in the 
detector, so both genuine high energy qq events, radiative returns and four- 
fermion hadronic events are selected for the hadronic lineshape. This gives the 
total sample: the background is subtracted on a MonteCarlo, using CC03, NC02 
and hadronic two-photon collisions. Clearly the above strategy is good enough 
for the present precision, but wrong in principle. Let us consider the relation 
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between (radiatively corrected) two-fermion (2F) and four-fermion (4F) final 
states in e + e~ annihilation at LEP 2. There are several components in the 
radiative corrections to fermion pair production: among them there is initial- 
state (or final-state) fermion-pair production. For definiteness consider e + e~ — > 
bb with radiation of an e + e _ pair [ ]l3| . The background is represented by the 
full four-fermion process, the so-called NC48 process, which is built out of 48 
Feynman diagrams. For studies around the Z resonance the default p"T[ was 
to included pairs from initial state and a cut was selected so that M{bb) > 
0.25 s. At LEP2 energies or higher one needs a more precise separation between 
radiative corrections to 2F production and real 4F events |l2| . 

We will denote the evaluation of any one-loop corrected cross section, e.g. 
e + e~ — ► bb as a 2F-calculation. By 4F-calculation we mean a tree level evalua- 
tion, e.g. e + e~ — * bbe + e~. Note that the soft pairs, 7* — > e + e~ are divergent 
in the limit of zero e + e~ invariant mass and therefore any simulation of very 
soft pairs with massless 4F-calculations is bound to produce wrong results. 

But also a massive 4F-calculation is not enough, because if pairs are soft 
enough we must include virtual pairs as well, and all e + e~ pairs are allowed 
down to M(e + e~) = 2m e . 

Also soft+virtual initial/final pairs in a 2F-calculation are not enough be- 
cause no upper cut is imposed on M(e + e~), so that all pairs compatible with the 
request M(bb) > (some thresholds) are accepted. Thus there is more than S+V 
pairs, there are many topologies for hard pairs and some of them require a finite 
m e also for hard pairs. Indeed in NC48 there are multi-peripheral diagrams 
which diverge for m e — > 0. 

The evaluation of a(qqX) requires |l3) 

A Include virtual+soft (up to some invariant mass A) I/F state pairs with 
a complete 2F-calculation. 

B The contributions to e + e~ — > bbe + e~ not in [A] are then included with a 
constraint on M(bb), but with no further restriction on M(e + e~), with a 
complete (i.e. fully massive) 4F-calculation. 

C The contributions to e + e~ — > 66e + e~that are already contained in [A], 
are included with M(e + e~) > A. 

Step B requires evaluation of the following 4F-processes: Fully hadronic, 
with at least one invariant mass passing the cut, semi-leptonic, with M(qjq i ) 
passing the cut. 

To summarize we may say that around the Z resonance the rate for real and 
virtual radiation is known and included in the existing calculations. Both 
are enhanced by large logarithms but they cancel to a large extent, leading to 
a small contribution to the inclusive decay rates. 

The complete evaluation of a(qqX) would be relatively easy if we could 
separate sub-classes of diagrams, e.g. primary from secondary production. For 
that it is necessary that the interference between them be zero, or very small or, 
at least non-singular. In the limit of massless fermions, singularities will arise 
from: 
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1. s-channel 7(g) 

2. t-channcl 7 with outgoing e ± lost in the beam pipe. 

This consideration suggests the appropriate strategy: classes of diagrams 
showing a mass singularity must be included through some analytical calculation 
which also accounts for O (a 2 ) virtual radiation, all interferences exhibiting 
mass singularities belong to this category while the rest, including most of the 
interferences, are accounted for by some (numerical) massless 4F-calculation. 

The goal of this paper is to investigate in more detail the class 2) introduced 
above, of which CC20 is a prototype. 

3 Kinematics and structure functions 

The inclusion of QED initial state radiation in e + e~-annihilation is based on 
renormalization group ideas and on factorization of mass singularities. The 
corresponding cross section may be cast into the following form: 

a(s) = / dxi / dx 2 @ cut V] D f (ari, s)D f + (x 2 , s) a fl (xix 2 s), (2) 

J " J " /=e+,e- 

where the structure function D£(x, s) is the probability density to find a parton 
f with energy fraction x. The restriction on the region of integration, given by 
8 C ut, reflects the presence of kinematical cuts. 

Let p be the four-momentum of the incoming electron in the laboratory 
system, 

1 m 2 
p=-^(0,0,/3,l), /3* = 1-A (3) 

2 s 

The electron, before interacting, emits soft and collinear photons. Let k = 
ki + k 2 + ■ ■ ■ be the total four-momentum of the radiated photons. Thus 

k= ^(1- x) (0,0,1,1), (4) 

so that k 2 — 0, as requested by collinear, massless, photons. Usually one can 
work with the massless approximation for the electron taking part in the hard 
scattering, thus an on-shell (massless) electron can emit a bunch of massless, 
collinear, photons and remain on its (massless) mass shell. But the electron 
mass cannot be neglected in the hard CC20 scattering and, after radiation, the 
electron finds itself in a virtual state having four-momentum 

p = p-k=^Vs(O,O,0-l + x,x), (5) 

with x being the fraction of energy remaining after radiation. As a consequence, 
the electron is put off its mass shell, 

p 2 = — ml + i(l — /?)(! — x) s <~ — xm\ for m e — > 0. (6) 
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When considering the whole process we introduce p± for the incoming e in 
the laboratory system. Once radiation has been emitted the momenta will be 
denoted by p± with 

P± = 1^(0,0,t(P-1 + x±),x±). (7) 

The total four-momentum becomes 

P = p+ +p- = (0,0, x_ - x + ,x^ +x + ) , (8) 

with a corresponding invariant mass 

P 2 = -x+x- s = s (9) 

In the following we will be able to discuss the effects of a correct treatment of 
QED initial-state radiation (ISR) on the processes under consideration. 



4 Weizsacker- Williams approximation for CC20 

The strategy for the calculation of the CC20 process will be as follows. First, 
we split the 20 Feynman diagrams of the CC20 family into the four diagrams of 
Fig. ^, characterized by the presence of a i-channel photon, and the rest 

CC20 = CC20 7 + CC20 R . (10) 

Then we introduce 9 C , the angle separating the SACC20 from the LACC20 
regions. The total cross section will be computed as 

|CC20<(m e )| 2 + |CC20 > (0)| 2 + 2 [cC20<(0)j f CC20£(0) + |CC20^(0)| 2 (11) 
where CC20 > (CC20 < ) implies 9 > 9 C (9 < 9 C ) and the argument m e (0) implies 



a finite(zero) electron mass. For the first term in Eq.(|ll[) we need an analyti- 
cal calculation which keeps m e =/= while the remaining terms can be treated 
numerically with the approximation of m e — 0. The square of the four dia- 
grams of Fig. H will be computed within the improved Weizsacker- Williams - 
approximation (WW), provided that 9 C is not too large. This application of the 
WW-approximation is very similar to the one applied in |L5|. 



w 
h 




+ 




Figure 2: The CC20 7 family of diagrams. 



The advantage of using the improved WW-approximation is in the possibility 
of performing an analytical integration over the momentum transferred to the 
photon, allowing to obtain the exact logarithmic enhancement as well as the 
first, constant, correction to it. 

The kernel cross section for the process CC20 7 (m e ), 



e + (p+)e (p-)->e (q-)u e (q + )u(k)d[k), 

can be written as 

_ 9 2 4 N, 



1 



(12) 



(13) 



where Q = p- — q- and N c = 1 for a fully leptonic final state and 3 otherwise. 
Furthermore sg is the sine of the weak mixing angle. In this equation, g?<I>3 is 
the phase space integral for the v e ud system, S + (p 2 + m 2 ) = 6{E)5(p 2 + m 2 ) 
and L, W are the leptonic tensor and the SACC20 tensor. A straightforward 
calculation gives 



L, 



l r 
2 



Q 2 - (1 - X-) m 2 e 6^ +p- fl q-„ +p- ly q- fi . 



(14) 



The four diagrams of Fig. g form a U(l) gauge-invariant set and therefore 
Q[iW^ v = QvW^j = 0. The VF-tensor admits a decomposition into three form 
factors 



/ 



P+ ■ Q 



w 2 



Q 2 



Q, 



Q a fP 

p+.q) p+-Q 



(15) 



Note that the W3 form factor gives zero contribution in this case. Actually 
the gauge invariance of CC20 7 poses a problem with a well-known solution: a 
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complete treatment would require the application of the Fermion-Loop scheme, 
but for our purposes it is enough to introduce a fixed width for the W, both for 
the s-channel and the i-channel, i.e. the fixed-width scheme. For a complete 
discussion we refer to the work in || . 

In the limit Q 2 — > 0, W^ v must be an analytical function of Q 2 . By requiring 
that Q 2 J d<£ 3 W^ = for Q 2 = one obtains 

W 2 (Q 2 ,y) ^W 1 (0,y)+o{Q 2 ) , (16) 
where we have introduced the variable 

y=^$-. (17) 
p+ -p- 

The WW-approximation is defined by the following equation: 

J rf$ 3 W^=>-\ Wi (0, y) / 7 (Q\ y, x 

Wi(0,y) = -5 J d<s> 3 w^ 2=Q (18) 

Wi is therefore proportional to the cross section for e + 7 — * v e ud (with real 
7) and / 7 is the photon density, which in the presence of QED initial state 
radiation reads as follows: 

h = -m\ (l + x _+2 l —^j+Q 2 (l-| + A). (19) 

Thanks to Eq. ( |l8| ) the Q 2 integration can be performed analytically. Note that 
the integrand is the sum of two terms, proportional to 

A (20) 
Q 2 Q 4 

To integrate over Q 2 we need the kinematics of the process which is specified, 
in the laboratory system, by 

p_ = ~sA(0,0,/3,l), q- = E f 08/ sin0, 0,j3 f cos 6,1), (21) 

with f3f = 1 — m 2 /E 2 and /3 defined in Eq.(||). Let Q(Q) be the momentum 
transfer with (without) inclusion of ISR, then 



Q 2 = x-Q 2 - (l-x-)m 2 + (l-{3)(l-x-)E f l3 f vGcos8 

= x^Q 2 — (1 — X-) (l — 2 —p cos# J m 2 + O (m 4 /s) . (22) 



,/i 

If we introduce the variable y, equivalent to the fraction of the electron energy 
carried by the photon in absence of ISR and defined by 

y = , (23) 

p+ -p- 



then the following relations hold for m e = 0, 



(l-cos^B/Vs, 2/= 1- (l + cos0) (24) 

V s 



Using this result in Eq.(|22j) one derives 



Q 2 = x_Q 2 - (1 - (y + ^ m 2 + O . (25) 

With ISR one uses y defined in Eq.(|l7]) and the relation between y and y is 
obtained from 

( Q 2 

2p + ■ Q = —s — 2x+p+ ■ q- + (x+ + x-)m 2 + (1 - x+) [y — H 

2 P+ -q- = (l-y)(2ml-s). (26) 
For finite electron mass the relations linking Q 2 ,y to Q 2 ,y read as follows: 

Q 2 

Q 2 = aQ 2 + bm 2 e y : y = c—+dy + e, (27) 

s 

where the a, . . . , e coefficients are 

2 

a = — (1 — x_) — -, b = —(l — X-), 
s 

1 — x+ m 2 , — £+ 1 + x + m 2 

c = — -, d = 



D s 1 D D 

1 — X- ( rn 2 \ „ . . m 



2 



n ... , D = -ar+z- + (X+ + z_) (28) 

V \ s J s 

We have introduced Q 2 and y because they are natural variables for describing 
the outgoing electron in absence of ISR. Indeed, one can show that 
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d 3 n 1 4- B : 

— p = 2Tvf3 f EfdE f dcos6 = tt — r -%- dQ 2 dy. (29) 
£/ 2/3 

The transition to hatted variables, to be used with ISR, is completed by deriving 
the jacobian of the transformation, 

dQ 2 dy =1 + (^pj dQ 2 dy. (30) 

Having specified the relevant variables we now proceed to deriving the bound- 
aries of the phase space. First we derive the boundaries for Q 2 . We start from 
the relations 

Q 2 = -2 m 2 e + (1-/3/3/ cos 0) E f y/s, 
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and introduce a new variable x denned by 



E f 



X 2 + m 2 e 



0f = 



2 2 

X -rnj 
X 2 + m 2 e ' 



From Eq.([jl|) a solution for x is 
11 + /J 2 



2 1 + fic 



i-v + J(i-y) 5 



1 — f3 2 c 2 m\ 
(1 + /3 2 ) 2 T 



(32) 



(33) 



where c = cos#. At zero scattering angle for the outgoing electron, c = 1, one 
obtains 



M ' 2 1 + 



1-y + O 



1 - y + J (1 - y) 2 - 4 



1 — (3 2 m 2 



{1 + I3 2 ) 2 s 



Inserting i?/ from Eq.(j32|) into Eq.(pi|) one derives 



Q 2 



9 1 JS 

-2ml + - — 
e 2 X 



(1 - 0c) X 2 + (1 + 0c) ml 



(34) 



(35) 



and, therefore the lower limit for the square of the momentum transfer is set by 

1 Js 



Q 2 (c=l) = -2ml + -V- (l-(3) X 2 + (l + f3)ml 

2 x L 



rn 



2 V 



o 



'"< 

s 



(36) 



so < Q 2 <Q 

,,2 



If we now require that 9 < 9 C , with 9 C <C 1, the limits for Q 2 are as follows 
-,2 ^ r& ^ n2 

where Xc is 

,2 /J2 



Xc 



to; 
s 



1 



+ -f (1 - y) + 



*e '"e />2 



giving an upper limit of integration for Q 2 

Q^ = Q^ + 7(l-y)^ s + Ol^,TO^ 



(38) 



(39) 
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The limits for Q 2 can be immediately translated into limits for Q, and one 
finds 



Q z < Q < Qi, 

Q 2 = m\ (x^y+l-x.) j^ + C[ — 



™4 



Ql = Q^ + \ X l{l-y)els + O[^,ur0;). (40) 

The photon flux-function, an essential ingredient of the WW-approximation, is 
now defined by 

^=[ Qc dQ 2 ±L. (41) 



Next we discuss the limits of integration for y and assume that the fermions in 
the final state are massless, apart from the electron. As we will see this can be 
the origin of new mass singularities. From this point of view ISR is inessential. 
Let us introduce variables p and k by 

l + (3 2 , , Q 2 + m 2 

such that the electron energy and scattering angle become 

E f = l(p + K)V^, pp f cos6=P—^. (43) 

Z p -T K 



For 8 = 0, after squaring the second relation in Eq.(|43|) and substituting Ef 
from the first one, one obtains 

K 2 -2±±^pK-p 2 +p 2 = 0. (44) 

In this way the allowed region of the phase space for the outgoing electron is 
completely specified. It is seen that for p > [3 one has k < or Q 2 < —m 2 , 
i.e. the square of the momentum transfer is not positive definite and crosses the 
zero independently of the finite electron mass. The variable k becomes negative 
for p > f3 or 

1 + /3 2 " s 2 

This simple fact is better illustrated by considering the process e + e~ — > e~X 
with X = {U e ud}. Let the cluster X be characterized by having four-momentum 
qx and mass M, i.e. q 2 x = —M 2 . The 2^2 process p+ + p- — > g_ + qx is 
described in terms of Mandelstam invariants 

s = -(p++P-)\ t=-(p_-«_) 2 , (46) 



y< K -W^r-^- (45) 
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so that Q 2 = —t and 



O 2 - m 2 + M 2 

V = ; + 2 ■ (47) 

s — 2 mi 



The physical portion of the phase space must satisfy the condition X > with 

X = i A (s, m 2 e , ml) A (s, m 2 , M 2 ) - s 2 t - m 2 e - M 2 + i (s + A/ 2 - to 2 ) 

(48) 

where A is the usual Kallen-function. When M = the condition X > is 
equivalent to 

t- <t< t+, (49) 

where one easily finds that 



27 to 6 

t_ — — ^, t + ~ -s, for m e — * 0. (50) 



Therefore, for M = 0, t is not negative definite. The amplitude squared is 
proportional to 1/Q 2 or to to 2 /Q and massless quarks induce a singularity, 
even for finite m e , if a cut is not imposed on the invariant mass M(«cZ)n. The 
singularity is, in any case, avoided by requiring a cut such that 



An upper limit on y is derived by considering again 

P+ + P- -* q- + qx , q x = q++k + k=p+ + Q. (52) 
Next we introduce the invariant mass of the quark-antiquark system, 

M 2 = -{p+ + Q f = ml -Q 2 + {s - 2m 2 )y, (53) 
and require the constraint 

Vs > m e + M, (54) 

equivalent to 

(s - 2 m 2 ) ?/ < s - 2m e Vs + Q 2 . (55) 
This inequality is satisfied for 

y<i-^, (56) 

following from the relation giving j/ ma x in terms of Q^n, 

(s - 2 to 2 ) y max = s - 2 TOeVs + Q 2 ^, Q 2 lin = Qo- ( 57 ) 

2 This fact was firstly pointed out in a private communication of A. Ballestrero. 
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Here Qq is taken from Eq.(|37]). The equivalent bound for y follows as 

0<1-^|, (58) 

With a cut on M{du) the singularity at Q 2 = (or Q 2 = with ISR) is 
avoided but we still have additional singularities. There are two multi-peripheral 
diagrams contributing to the CC20 process e + e~ — ► e _ 77 e /i/ 2 , the last two in 
Fig. H When Q 2 = 0, i.e. the electron is lost in the beam pipe, and the 
(massless) /i(/2)-fermion is emitted parallel to the (quasi-real) photon then the 
internal fermion propagator will produce an enhancement in the cross section. 
Taking into account a him 2 from the photon flux-function, 3 options follow: 

1. to consider massive fermions, giving a result proportional to mm 2 lnm 2 , 

2. to use massless fermions, giving instead In 2 m 2 , 

3. to introduce an angular cut on the outgoing j\ and f 2 fermions with 
respect to the beam axis, 0(/i,/2) > C ut, giving In m 2 In 6 C ut ■ 

The first option is clean but ambiguous when the final state fermions are 
light quarks, what to use for m u , mdl The second one presents no problems for 
a fully leptonic CC20 final state but completely fails to describe quarks, as we 
will show in discussing QCD corrections. The last option is also theoretically 
clean and can be used to give differential distributions for the final state jets. It 
is, however, disliked by the experimentalists when computing the total sample 
of events: hadronized jets are seen and not isolated quarks. Even if the quark is 
parallel to the beam axis the jet could be broad enough and the event selected. 

These events are also interesting since they correspond to a situation where 
the electron and one of the quarks are lost in the beam pipe, while the other 
quark is recoiling against the neutrino, i.e. one has a totally imbalanced mono- 
jet structure, background to new particle searches. 

The singularity induced by massless quarks in e + e~ — ► e~T> e ud can only 
be treated within the context of QCD final state corrections and of the photon 
hadronic structure function (PHSF) scenario. We will come back to the problem 
later in the paper. In the next section we discuss, instead the cross section for 
the sub-process e + 7 — > V e ud in the two regimes, massless and massive quarks. 
It will be seen that one can use different parametrizations for the corresponding 
phase space, depending on the presence of kinematical cuts. 

5 The sub-process e + 7 — ► V e ud with cuts 

As a consequence of Eq. (|l8|) , the WW-approximation, we write the result for a 
CC20 cross section as the convolution of the photon flux- function of Eq. (|d[) with 
the cross section for e + 7 — * v e ud. The process e + (p+)j(Q) — > V e (q+)u(k)d(k) 
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is illustrated in Fig. || and is described by the following invariants: 



P+ ■ Q = 




1 

— ys, 
2 y 




1 

q+= 2 








Q -q+ = 


1 
2 


K -i 




1 

= -t, 
2 ' 


Q 






q + -k = 


1 
2 


C-, 












k-k = 














(59) 



The process takes place at energy ^fys and there are only five linear-independent 
invariants, which we choose to be ys and r, x\,x 2 and z defined by 

t~rys, s' — ^—xiys, s — x 2 y s, k~ = zy s. (60) 

The phase space can be computed in terms of the following object: 

= fs 2 J d 4 kd 4 kd 4 q + 5+(k 2 )5 + (k 2 )5 + (q 2 + )5 4 (p + + Q-k-k-q + ) 
x S (2 k ■ k + x 2 y s) S (2 k ■ k + 2q + ■ k + x 2 y s) . (61) 




The integration is most conveniently performed in the system where 

P = P+ + Q= (0, 0,0, VP) . (62) 

Moreover, let k be along the positive z axis with Q in the x — y plane (and polar 
angle denoted by ff) and let q+ be described by angles ip, <fr. Then we find, 

J, 



dx\dx2 



y 2 s 2 e( Xl ) d 4 kS+(k 2 )s(2E d ^y~§~x 1 ys) 



J = 6(1 — x 2 ) — / dcostp 

8 J -i 

x 6 (-xi(l-x 2 )(cos<i/)- l)ys + (xt -x 2 )n - : ) ■ 
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It is more convenient to introduce 

d 3 <f> 3 

dx\dx2dz 

J 



y 3 s 3 v{xi 



i) / d k5 + (k 2 ) 5 \2Edyys — xiys) J, 



;i-s2) 



1 — X2 



+ 1 p2ir 

d cos V> / d</) 
-l Jo 



x (5 (((1 — X2) sin sin ip cos + (1 — £2) (cos 9 cos V> — 1) — 2 z) y s) 

(64) 



x * ( 2 ~ ^2)(cos^ - l)ys + (ati - 



If we now take into account that 



cos 9 = 1 + 2 — , 



the final result follows 



2 * - 

7T 7/ S 

ii$3 = Qdx\dx2dzdr, 

Ax\R 



(65) 



(66) 



9 = {JJ 0(3^)0(1 -3Gi)}0(a;i -x 2 )6>(z+ - z)0{z - Z-)6{-t)Q{t + x x ), 



where gives the boundaries of the phase space and where we have introduced 



R z 

cos tp 



-4z 2 + 4 (1 - a; 2 )(cos0cos^ - 1) z - (1 - x 2 ) 2 (cos0 - cost/)) 2 , 



2 X2 — x\ — 



X\{\ - x 2 ) 
and where z± are the roots of R 2 = 0, i.e 

1 — X-2 



z± 



(cos 9 cos ip ± I sin sin?/' I — 1) . 



(67) 



(68) 



The cross section will be computed with a cut on the invariant mass of the 
ud-pair, i.e. M 2 (ud) > sq giving 



so 

— < x-x < 1, ys > so- 
ys 

Starting from the four diagrams of Fig. we derive 
W = 



(69) 



9 s e 



f Wu + I A t | 2 W^ 22 + I A, I 2 W 33 + y|A s | 2 W M 



f (Wu + W 13 + Wu) + y (w 23 + W 2i ) + y |A S | 2 W< 



:S1 



(70) 
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where the propagators are defined, in the fixed width scheme, by 
A s = -x 2 y + n 2 w - i 7 w Mw A * = ~( l ~ x z + z )v + A w ~ *7wMw> ( 71 ) 

with T w = 7 w v / ^ and fj, w — M w y/s. The Wij represent the i — j interference 
of diagrams in Fig. || and are given in the following list: 

Wn = t(1 + zx\ + 2z - 2xix 2 + xi + z 2 ) 

+ t 2 (1 + z — x 2 ) + zx\ + zx\ + 2x\X2 + xiz 2 — 2x 2 x\ — x 2 + x\, 



W22 = 

= 

W 12 = 



yr(l - xi) 
1 r 



--(l + Xl ) + -r 
9v 9 



--(l-aji , 
9 T 



9 v 



1 Xi 

9~v 



ReA, 



r(— 2 + z:Ei — 2z + 2xix 2 ) 



+ t 2 (-1 + x 2 ) - zx\ - 2xix 2 + 2x 2 x\ + x 2 



W 
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^ReA f 



— (zx\ — 2zxf — 2xix 2 + x\ + 2x2X\ + x 2 — 2xf) 

T 



2 , .2 2 4 2 „ 



1 rr 1 

W^i4 = -ReA s -(-1- 2x1x2 + 3xi -2a;?) + -r(l + 2z) 
3 L<y 3 

1 „ s 1 1 r n n 2 s 2 

+ 77 — (1 -x x -x 2 ) + --{2xix 2 - 2x 2 x x - x 2 ) + -zxi 
3 v av 3 



W; 



2:s 



= — ReA s ReA f 
3 



•^1 „ _ N 2 , 4 4 2 



— +-xi--x 1 
t 3 3 



4 7^M^ /, „ ^2 



ia;i(l-si)^ 



+ 37^^1(1-^1) 



1^24 



1 r 7- 1 l T 2 

= -ReA s ReA f -(1 + 2xhe 2 + xi - 2x\) + -r(-2 + n) + - — 
3 Lu 3 3 u 

x (1 - xi + x 2 ) + --(—2xix 2 + 2x 2 Xi + x 2 + 2x\ - 2x\) - \x x 
6v 3 



+ zlivi-i 1 + 2x x x 2 +xi- 2x1) + T } li^ w T-{-2 + xi) + ^lfi 2 w 
3 v 3 3 v 



X (1 - Si + X 2 ) + 

3 v 



(-2x1X2 + 2x 2 x\ + x 2 + 2x\ - 2xJ) - -"fl^Zxi, 



2 1 2 1 

W34 = - — (-2xix 2 + 2x 2 x 2 +x 2 ) + --xi(l-2xi) 

9 TV 9 T 

2 t 2 1 2 

+ --(-1 + Xi + x 2 ) + --(1 + 2xix 2 - 3xi + 2x\) - -xi 
9 v 9 v 9 



(72) 
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where we have introduced a new variable v, 



V = 1 + T + Z. 



(73) 



There remains the problem of angular cuts for the outgoing quark (antiquark). 
These cuts should be applied with respect to the beam direction and should be 
expressed in terms of the variables describing the sub-process. For CC20 7 and 
for 9 C <C 1 the e~ is lost in a narrow cone around the beam, so that we can use 
the approximation Q = yp- . In the laboratory system we have 



p± = ~^(o,o,tM). 

Let E u be the energy of the outgoing u-quark, so that 

k = E u (sin6> u ,0, cos0 u , 1) , 



(74) 



(75) 



where 9 U is the scattering angle of the u-quark with respect to the incoming 
electron. One finds 



COS* 



P+ 



p+ • k — p- ■ k 
p + ■ k + p- ■ k 
P+ ■ k u' 
x + 2x+' 



V- 



k 



t 



Q-k = 

X-y 2x-y 



(76) 



The condition | cos8 u \ < C becomes, in terms of invariants, 

{l + C)x-yx!+ (1 + C)x m y+(1-C)x+ t < 0, 
(1 - C) x-yx! + \{l - C) x m y + (1 + C) x+] t > 0. 



(77) 



With k — Ed (sm9d,0,cos9d, 1) we derive two additional conditions similar to 
those of Eq. (|77j) but with v! — > t',t — > u and reflecting the cut | cos 9d\ < C. 

Eq.([7C|), in conjunction with Eqs.(]T3|, p^[ ) and Eq.(|67|), allows us to compute 
the cross section within the WW-approximation. 



6 The fully extrapolated sub-process e + 7 — > v e ud 

Our goal is to compute the CC20 cross section without any kinematical cut, 
apart from imposing that M(ud) > y/so. We have seen in sect. ^ that there 
is a mass singularity in the total cross section when the quarks are assumed to 
be massless. Here we compute again the cross section with finite quark masses. 
Let us consider again the process e + (p + )j(Q-) — > v e {q + )u(k)d(k). Moreover 
let Q+ be defined as Q+ = p+ — q+ and let Q be Q+ + Q- . In order to compute 
as required by Eq.([l^), we neglect for the moment ISR and introduce 
three master scalar-integrals: 

/„= f d 4 k d 4 fc^+(fc 2 + ml) S + (k 2 + m 2 d ) S 4 (Q -k-k) 1 » , (78) 
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with n = 0, \2. They are easily evaluated in the system where Q = (0,0, 0,/z) 
and where the photon four-momentum is Q_ = X (0,0, 1, 1). One immediately 
finds 

X = 9±^. (79) 



If we start with I 2 then we obtain 

S+(2uE„ - a 2 -to 2 +TO 2 /! 

(Q_ • k) 



h = I d 4 k 5+ (k 2 + ml) 8+ (2 iiE u - M 2 - to 2 + to 2 ) 1 — 
J (Q- ■ 



for m u -> 0. (80) 



4^X 

7T £ 7T Q 2 1 

4 M^ 2 m 2 ~ ~8 (Q + .Q_) 2 m 



Similarly we evaluate I\ as follows: 

Ji = J d 4 k8+(k 2 +m 2 u )6+(2^E u -^ 2 -m 2 u m 2 ) 



X (kcos9- E u ) 



7T TO„ 7T TO 2 

= m - ~ In — £r tor m u — > 081) 

2 MX £?„ + ^/^ _ m 2 4Q+-Q_ -Q2 

Finally, for 7 one gets 

= f • (82) 

The complete result for the cross section follows from squaring the matrix ele- 
ment, 

W^ = ^Kt(-it + m u )K li (it + m d ) (83) 
The function 1Z is 

K» = 1 a 1+ nl a P(s)P(t)+ri+^aP(s)P(e) 
+ l^Q-fi-im^ 7 "7+ K P(t)P(f) 
+ l a l+ + i m d ) 7 X P(t)P(f)- (84) 

Furthermore we have 7+ = 1 + 75 and 

K\ a = V f ,p a v(p + )~f p ~f + v(q + ), 

R>1 = -Quv(P+hv.l+v(q+), 

K-l = Qdv{p+)lnl+v{q + ), (85) 
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where V^p is the corresponding triple gauge-boson vertex, Q u (Qd) is the up- 
(down-)fermion charge and the propagators appearing in Eq.(84) are 



p(a) 



P(e) 



1 



(k + kY + Ml-iT w M w 
1 



P(t) 



(p+-Q-) + m l 

p(/) = — t> n/o = 

(Q_ - fc) + m£ 



(Q_ -k) + m 2 



(86) 



Next, we have to integrate over the phase space. The integration over fc, k can 
be performed by introducing additional integrals: 

= J d 4 kd 4 kS+{k 2 + ml)S+(k 2 + m 2 d )S 4 (Q-k~k) 
k^ 1 ... k^ 



(Q_ • fc) 

d 4 fc d 4 k5 + {k 2 + ml) 5+{k 2 + ml) 6 4 (Q -k-k) 
k^ 1 ... k^ 



(Q--k) 



(87) 



All these integrals can be reduced to the scalar form factors. Quarks masses are 
kept only in front of I2 and the reduction gives 



3n 



"(?) = 



Q l "Q v Q a hnMd) + Q-Q-Q-hnMi) + (Q-Q"Q a + Q»Q-Q a 
Qi*Q"Q1)hnM + (Q"Q-Q- + QtQ v Q- + Q^Q-Q^huMq) 



1 30 



h 

4 



^Q v Q a - {Q' 1 8 va + Q v 8' ia + Q a 5> lv ) 



6 



= Q^hnM + Q-Q-hnMi) + (Q' l Q- + Q"Q-)hnMi) 



1 20 



h 

3 



h 
2 



(88) 



where q — u,d and where all the form factors can be reduced to linear combina- 
tions of the master scalar integrals of Eq.((7l). After a straightforward algebra 
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one obtains 

Ji2,n(g) = 

J12, 12(g) = 

hi, 11(g) = 

hi, 12(g) = 

J22, 24(g) = 

-^22, 22(g) = 

^21, 24(g) = 

J21, 21(g) = 

^21,23(g) = 

^21,22 (g) = 

-^32,36(g) = 

-^32, 32(g) = 

^31,36(g) = 

J31, 35(g) = 

^31,31 (g) = 

hi, 33(g) = 

^31,34(g) = 

^31,32 (g) = 



0, 

h(q)[l 
h 

Q+-Q- 

h(q)-2 



1 



2 Q+ • Q. 



Jo 



0+-Q- 



1 0+ 



2Q+ 0- 



0, /22.21(g) = 0, /22.23(g) = 0, 



h(q) 
1 



-/r 



2 

1 /o 



1 + 
1 



1 

2Q+0_ 

1 2 + 



2Q+ -Q_' 

1 Jo 

2Q+ • Q_ 

/i(«)-3 



2Q+ O- 
1 o 2 + 



1 + 



/o 



2Q+ Q- 



0+-Q- 



i + i Q 



20+ -0- 



0, /32.35(g) = 0, /32.31(g) = 0, /32.33(g) = 0, /32.34(g) = 0, 



1 T 



1 



0+ 



(0+-0-) 
1 Qi ^ 



3 

+ - 



0+ 



4(0+-0-) 2Q+-0- 



20+ -0- 



1+1 



?+ 



20+ -0- 



3 0+0-' 
1 Jo 

6 0+ -0- 

1 Jo 



1 



ft 



2Q+ 0- 



3 0+0- 

Ji(g) 



1 + - 



1 0+ 



i + i- 



2 0+0- 

0+ 



Ql 



3 



- 11 Jo 



(0+-0-) 3 4 (0+-0-) 2 ' 2 0+-0- 
1 Qi 1 01 , 1 Ql 



24 (o+-o-r 4 (o+-o-r 2 (o+-o-r 



+ 



1 1 

3 0+ -0- 



(89) 
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ISR is restored by changing in the previous equations Q, Q± into Q, Q±. The 
kernel cross section for the process, to be convoluted with the e^- structure 
functions, is therefore written as 



dQ 2 _dy-L / d^W, 



._ £4 * 

(2tt) 8 64 s J ' ' J ~ Q2 =0 

J d$ 3 W^ 2=o = ^N c ,J dQ 2 + dy + \A\ 2 , (90) 

where A is the amplitude for the process, function of Q\, y±, with 

Q+=p+-q+, y+ = -2 — ;— , y- = - — — • (91) 

y-s p+-p- 

Starting from the original d$3 we have been able to perform the k, k integrations, 
with the help of Eqs.(^7|-^9|), arriving at a twofold, dQ 2 dy + , integral. For the 
purpose of integration it is more useful to change variable from Q 2 ^ to x, defined 

by 

Q\ = (y+-x)y.s. (92) 
The limits of integration and the jacobian of the transformation are: 

t^-x < x < y + , ~~~~ < y+ < 1, dQ 2 + dij + = y-Sdy + dx. (93) 

y-s ys 

Before giving the complete expression for \A\ 2 and computing the cross sections 
we have to answer the question of what to do with the light quark masses. The 
following section is devoted to a clarification of the origin of these additional 
mass singularities. 



7 QCD corrections 

We have already indicated that, for massless quarks, the cross section for e + e~ — * 
e~v e ud is dominated by two large logarithms. One originates in the limit of 
small scattering angle of the outgoing electron. The other comes from the prop- 
agator of the internal light quark in the multi-peripheral diagrams. Another 
way of looking at it is to reconsider the integral I\ of Eq. (|8l|) and to evaluate it 
for m u — wid — and arbitrary Q\. One obtains 

^ 0) = [ d 4 kd 4 k8 + {k 2 )5+{k 2 )5\Q~k-k) 1 



(Q--kf 



= *Tn «+-V- — (94) 
4VA Q+-Q- + VA 

with A being a Gram's determinant, 

A = (Q+-Q_) 2 -Q\Q 2 _. (95) 
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Therefore, for small scattering angles of the outgoing electron, the integral be- 
haves like 



4 0) ~ A *„ ^ Q2+Q ~ a , for Q 2 _-0. (96) 

Note that (Q- — k) 2 appears in the internal quark propagator of the multi- 
peripheral diagrams. For values of Q 2 _ small enough, the lower limit for {Q- — 
k) 2 becomes much smaller than Aqcd, well beyond the limit of applicability of 
perturbative QCD. This fact has many similarities to the inelastic ep scattering, 
see the work in Q . 



Figure 4: Example of resolved type processes in CC20 



So far, QCD corrections to the CC20 process have been applied within the 
context of naive QCD 0] or with a complete O (a s ) calculation |l7j which 
assumes a point-like coupling of the photon to quarks. However, the large log- 
arithm of Eq.(p6[) receives contributions from any order in a s from multiple 
gluon radiation. The latter creates a series of extra quark propagators, each 
yielding an extra power of the logarithm compensating the additional power 
of a s . This fact is discussed in and in |fL8f . Therefore, logarithmically en- 
hanced terms, of order a™ In™ (m 2 / Q 2 ) appear at every order in the perturbative 
expansion and, since the logarithm is large, the perturbative series does not con- 
verge quickly. Fortunately, this difficulty can be obviated, at least in principle. 
A formalism exists to sum these logarithms to all orders in perturbation theory, 
see @. 

The point is that a photon has a point-like coupling to the quark-antiquark 
pair only for sufficiently high virtuality. On the contrary, for small electron 
scattering angle, the photon is quasi-real in our CC20 process and behaves like a 
meson. This fact and its consequences are well-known in other processes, like the 
total jp cross section. At low photon virtualities one also expects contributions 
from the partonic constituents of the photon. The two contributions form the 
so-called resolved photon component, see Fig. |4|, which is usually added to the 
direct one, computed to fixed order in perturbation theory, where the photon 
is treated as an elementary particle. The resolved photon component is given 
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in terms of the photon, hadronic, structure function. There is a well-known 
subtlety in adding direct and resolved components. The direct component, 
evaluated at some fixed order in a s from all corresponding Feynman diagrams, 
contains singular terms that are already re-summed in the resolved component. 
The correct result flfl is schematically represented in the following equation: 



u = I dy- T. 



«res i ~dir.sub 



(97) 



where the superscript sub indicates that one must subtract the terms responsible 
for the large logarithms in the direct photon component and where <r 7e is the 
direct or the resolved cross section for r ye + — > V e ud at fixed y_. In the resolved 
cross section the photon interacts via the quark or the gluon component in its 
structure function. Therefore one has Q 

V™ =J2 J (Jl, M) *ie-jet. (98) 

here F i7 is the PHSF and M is the factorization scale. As shown in Eq. (|97|) the 
cross section for the CC20 process, e + e _ — > e~z7 e jet follows from the cross sec- 
tion for e + 7 — > F e jet by applying the equivalent photon or Weizsacker- Williams 
approximation which factorizes the flux of quasi-real photons emitted by the e~ 
from the interaction rate between the positron and the photon assumed to be 
real. 

The introduction of a resolved component for the photon is a familiar topic 
in 7p or 77 scattering. Here the situation is slightly different. The small virtual- 
ly of the photon is only needed when CC20 is a background to the high-energy 
hadronic lincshape or to Higgs boson searches or for single W production. The 
cross section for CC20 7 is obtained starting from four Feynman diagrams, two 
being single W resonant and two being multi-peripheral. In single W produc- 
tion, where one applies a cut |cos#(e~)| > 0.997, the QCD corrections are 
important and it appears difficult to obtain a precise prediction for the total 
cross section without summing the large logarithms into the Fi 7 distribution 
function. 

For the hadronic lineshape, on the other end, hadronic events are selected 
based on final state particle multiplicity in the detector. This gives the total 
sample but more interesting is the high-energy M 2 {ud) > so sample. The rela- 
tive dominance of the multi-peripheral diagrams in CC20 7 is larger in the total 
sample but not necessarily in the high-energy one. Therefore the uncertainty 
associated with the use of the PHSF calculated at the zeroth order in a s , i.e. 
in the Born approximation, is less relevant if we apply a strong invariant mass 
cut. 

Our strategy, for the moment, will be to use the parton model result, i.e. 
zeroth order in a s , and to cure the ill-defined massless limit by replacing the 
quark masses with a factorization scale M , in our case m u — — > M . The 
total cross section for the CC20 process will, therefore, depend on the scale M. 
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The amplitude squared, to be inserted in Eq.(fOT), becomes 



4 y_ 



2„,2 ,,2 



3|A. 
1 V- 



l2 V Iw^ 



Y y% 



3 I A, I 2 yl 



2 y- 3 2 xl-y. 



3|A S 



2 yt 4 x 2 l-y 



3|A S | 2 F y3 



2 - 4y+ - $ + y 2 + 2x(l - y+)(-y+ + x) 
y+ + x(-3y+ + 2x) 
y\ + x{3y+ - 2.x) 



Y M 



1 V- T v 2 



-5f+ -f+ + 2y% + x(15f + - 6f + + 3 - 20xy+ + 8xy: 



9|A t |2 yj 

Axyl + 10x 2 - 4x 2 y+ + 2x 2 y\) 



4 y_ 



,,2„,2 ,,2 



x l - y+ 



3 I A*. 
1 y_ v 2 



2 V ^w^w v ~4 



y 2 + x(-3y+ + 2x) 
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9|A t |2y4 

32xy 2 - 8xf + + 30x 2 - 18x 2 y + + hx 2 y 2 + ) 



~y+ — 50x^4 



2 y 2 



,,3, ,2 



3|A t 



s i-i/-t 

y y3 



-y 2 + x(3y+ - 2s) 



2 y^ 4 i 1 - y+ 



-y3 + a;(4y^. - 5a$+ + 2x 2 ) 



3|A t |2" y f + 
The propagators, in the fixed width scheme, become 

1 , 1 



where 



-xy-v + ^ w -ij w ^ 



v = x + x-, Y 



A, 



(x - y+) y-v + n 2 w — i j w fi l 



1 ,,2 ■ 

w 



L = bx- 



M 2 
xy-s' 



(99) 



(100) 



(101) 



and y±,x are given in Eqs.(^l| |92]), u w , j w after Eq.(fnj). With \A\ 2 at our dis- 
posal we can use Eq.([X)|) and Eq.(|93|), apply the convolution with QED structure 
functions, and derive the total cross section. 
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8 Numerical results and conclusions 



In this section we present all relevant numerical results for the CC20 processes 
as computed by the FORTRAN program WTO version 2.0 The chosen 

setup is specified by the following list: 



a/s = 186 GeV, M w = 80.39 GeV, M z = 91.1867 GeV 



(102) 



Naive QCD is not introduced which implies, in particular, that the W width 
is included without QCD corrections. For our setup this results into T w — 
2.0459 GeV. The QED radiation is included by means of the structure function 
approach (in the so-called /3-scheme ]2(| ^l], [2^] ) . First we consider a cut on the 
scattering angle of the outgoing quarks with respect to the beam axis, 10° < 
8 q < 170°. We also fix a lower cut on the invariant mass of the ud system, 
M 2 (ud) > 0.01 s. According to the procedure described in Eq. (|i"l|) we introduce 
a separating angle 9 C and compute the following cross sections: 

cr< or |CC20<(TO e )| 2 for 9 < 6 C , 
ct> or |CC20>(0)| 2 for 9 C < 9 < it, 



o-<mt or 2 



CC20<(0) CC20^(0) + |CC20^(0)| 2 for for 9 < 9 C 



Our reference values will be 9 C = 0.3°, 0.4° and 0.5°. We have verified that 
c<int is completely negligible for our choice of the separator 9 C so that the total 
is safely given by the sum a < + cr> . Indeed we find cr<i n t = 6 -j- 5 -j- 3 x 10~ 5 pb 
for 9 C = 0.5° ^ 0.4° ^ 0.3°. Always for 10° < 6 q < 170° we find for a < /a > the 
results shown in Tab.(nl). 



0c [Deg] 


^< 


^> 


(T< + CT> 


0.3° 


0.0527 


0.6316(9) 


0.6843(9) 


0.4° 


0.0554 


0.6289(7) 


0.6843(7) 


0.5° 


0.0575 


0.6269(6) 


0.6844(6) 



Table 1: Cross section in pb for the process e + e — > e v e ud, for 10° < 9 q < 
170°, q — u,d as a function of 9 C . 

Tab.(Q) clearly shows that there is a smooth matching of the two components, 
< (m e ) and > (0) at 9 = 9 C . This result justifies the application of the WW- 
approximation in the narrow cone around the electron axis and we conclude by 
quoting the following result: 

0-CC20 (186 GeV, 10° < 9 q < 170°, M 2 {ud) > 0.01 s) = 0.6843(4) pb. (103) 

We have also varied the angular cut on the outgoing quarks, keeping M 2 (ud) > 
0.01 s and 9 C — 0.5°. The latter is fully justified by the tiny dependence of the 
cross section on the separating angle 9 C . The results are illustrated in Tab.(0). 
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O q [Deg] 


<?< 


^> 




5° 


0.0607 


0.6404(6) 


0.7011(6) 


6° 


0.0601 


0.6384(6) 


0.6985(6) 


8° 


0.0588 


0.6333(6) 


0.6921(6) 


10° 


0.0575 


0.6269(6) 


0.6844(6) 



Table 2: Cross section in pb for the process e + e — > e v e ud, for 6 C = 0.5° as a 
function of 9 q . 

Next, we consider the total CC20 cross section, without angular cuts on the 
outgoing quarks and with massless quarks. As explained in Sect. [7], where 
we have discussed QCD corrections, the resulting cross section depends on a 
factorization scale M. For M — 1 GeV we find the results of Tab.(^). From 



0c [Deg] 


o< 


^> 




0.3° 


0.0850 


0.6503(8) 


0.7356(8) 


0.4° 


0.0881 


0.6473(7) 


0.7354(7) 


0.5° 


0.0905 


0.6451(6) 


0.7354(6) 



Table 3: Cross section in pb for the process e + e — ► e v e ud, for the factoriza- 
tion scale M = 1 GeV, as a function of 9 C . 

Tab.(||) we derive the CC20 cross section in a fully extrapolated setup. 

0-CC20 (186 GeV, M = 1 GeV) = 0.7354(4) pb. (104) 

We have also investigated the dependence of the cross section on the factoriza- 
tion scale M. With 6 C = 0.5° and M 2 (ud) > 0.01s the results are presented 
in Tab.fl), showing a mild dependence on M of the total. Finally we have 



M [GeV] 




ct> 


cr< + (T> 


0.1 


0.0908 


0.6451(6) 


0.7359(6) 


1 


0.0905 


0.6451(6) 


0.7356(6) 


10 


0.0902 


0.6451(6) 


0.7353(6) 


100 


0.0900 


0.6451(6) 


0.7351(6) 



Table 4: Gross section in pb for the process e + e — > e v e ud, for 6 C = 0.5° as a 
function of the factorization scale M . 

analyzed single W production with 9(e~) < 0.5° in WW-approximation. In 
Fig. H we have reported the M(ud) distribution for 10° < 6 q < 170°. In order 
to understand the role of the different components we have plotted the distri- 
bution with and without the multi-peripheral component. It follows that this 
component dominates at low invariant masses, while above M(ud) « 70 GeV it 
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Figure 5: Invariant mass distribution for e + e — > v e ud showing the effect of the 
multi-peripheral component . 

is practically without influence. A final comment is devoted to the validity of 
the WW-approximation. Note that we only use this approximation in a narrow 
cone around the electron axis, typically 9 < 0.5° = 8.73 mrad, and the complete 
calculation outside the cone. Corrections to Eq.(|l8|) of O (Q 2 ) have been dis- 
cussed in [ pl| where it has been shown that, after integration, the cross section 
receives additional contributions proportional to powers of Q 2 / (2 Q ■ Since 

Q - <-^e (los) 



2Q-p+ 4 s 

and, in our case, ^ = 186 GeV, s = 0.01s we find Q 2 /(2Q-p + ) < 1.9 x 10~ 3 . 
Therefore non-factorizable corrections are formally negligible. 

The mild dependence of the total cross section on the factorization scale 
M can be understood from Tab.(||). Here, for 6 < 0.5°, we have reported: 
1) the total cross section , i.e the sum of the single-resonant and of the multi- 
peripheral, non-resonant, components; 2) the multi-peripheral component alone. 
Therefore the non-resonant, M-dependent, terms are strongly suppressed and 
the M-dependence has little influence on the total, justifying our approximation 
of working at zeroth order in a s . 
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M [GeV] 


„TOT 

a < 


of 


0.1 


0.0908 


0.0009 


1 


0.0905 


0.0007 


10 


0.0902 


0.0005 


100 


0.0900 


0.0003 



Table 5: Total Cross section and Multi- Peripheral component in pb for e + e — > 
e~V e ud, for 6 C — 0.5° as a function of the factorization scale M . 

The experimental Collaborations at LEP are now recording and analyzing 
a sizeable fraction of events with four fermions in the final state. Outgoing 
electrons represent a notorious problem because of the presence of i-channcl 
photons interacting with W bosons or coupling to quark-antiquark pairs. 

The collinear limit forbids a calculation where the massless limit for fermions 
is assumed from the beginning and, in turn, this may induce numerical instabil- 
ities in computing the total cross section, even for a fully massive MonteCarlo. 
It should be mentioned also that the majority of the MonteCarlos used in the 
analysis are built in the massless approximation and because of that the total 
cross section is not available. 

We have suggested a simple but realistic solution based on the use of the 
Weizsacker- Williams approximation, to abe applied in a narrow cone around 
the beam axis. The large logarithms, ln(mg/s), are correctly described by our 
numerical solution. Furthermore the improved WW-approximation that we are 
usingis valid beyond the leading logarithmic approximation, as explained in 
Eq.(|l8|), and correctly integrates also the m^/Q 4 terms present in the photon 
flux. 

The correct treatment of the kinematics, accounting for the introduction of 
QED initial state radiation, is also emphasized. We have derived a version of the 
flux-function which describes quasi-real photons emitted by the electron after 
QED radiation. A second logarithmic enhancement in the cross section, arising 
from internal fermion propagators, is also described and a link is established 
with the familiar examples of "/p or 77 scattering. Finally, several numerical 
results are shown, proving the goodness of the adopted solution. 
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